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Abstract. In this paper the inverse scattering problem for the nonstationary 
Dirac-type system on the whole plane was considered. A nonlinear evolution 
sytem of equation related to nonstationary Dirac-type system is introduced 
and the solviblity of this sytem using the 1ST method is studied. 



1. Introduction 

As is known, every linear scattering problem determines a class of nonlinear 
evolution equations for which the inverse scattering transform (1ST) method is 
suitable to integrate these equations. The investigation of the nonlinear evolution 
equations by the 1ST method has been stated in [1, 2]. The essence of this method 
is to represent this nonlinear equation as a Lax equation 

AT 

(0.1) ^ = [L,A] 

by using the Lax pair L and A. The equation in the form of (0.1) allows to 
investigate evolution of scattering data instead of the evolution of the operator L, 
if the inverse scattering problem is investigated for the operator L. 

The 1ST method for the nonlinear evolution equations with 1+1 dimensions 
(one space and one time dimensions) have been described in a variety reviews and 
monographs. The generalization of the 1ST method to nonlinear evolution equations 
with 2+1 dimensions (two space and one time dimension ) has been developed in 
monograph [3]. In this direction one can also notice the papers [4 - 7], in which 
the inverse scattering problem for the first order systems is studied and this inverse 
problem was applied to integration of N- wave interactions, Davey-Stewartson and 
Kodomtsev-Petviashvili equations . 

Let L = — M, where M is an ordinary differential operator in x. Then 
equation (0.1) is rewritten as (see [2]) 



(0.2) 



— -M — 

dy ' dt 
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and it is a relation for ordinary differential operators. Some nonlinear evolution 
equations in 2+1 dimension arc presented in [8] by using commutativity condition 
(0.2) when M is an scalar coefficients ordinary differential operator in x. 

This article is organized as follows: In Section 2, we consider the inverse scat- 
tering problems for the L = ^ — M, with M = a-^ + Q(x, y) in the whole plane. 
We introduce the scattering data for the operator L in the plane (the minimal 
information for a unique restoration of matrix-function Q(x,y)). All the results 
in Section 2 are the matrix generalization of the L. P. Nizhnik's results for two 
component nonstationary Dirac system on the whole plane in [9] (see [4] also). 
Taking into account that this generalization is not difficult, so we do not give the 
proofs. In Section 3, we introduce some system of nonlinear evolution equations 
with 2+1 dimensions related to nonstationary Dirac- type systems and we study the 
solvability of this system using the 1ST method. In [10], the 1ST method is applied 
to integrate the two component nonlinear Schrodinger equation, using the inverse 
scattering problem for stationary Dirac-type system in whole line. 

2. Inverse scattering problem for the Dirac-type system in the plane 

Let us consider a system of first-order partial differential equations (PDEs) 
dtp dtp 

(1.1) L(tp) = — -a—-Q{x,y)tp = 0, - oo < x, y < +oo, 

where 



1 qi 

1 ,Q = q 2 
-1 J [ q 3 q A 

with the measurable complex- valued functions qi (i = 1, . . . , 4) . 

We give fundamental results on the inverse scattering problem for the system 
(1.1) in the case where the coefficients qi( i — 1, ... ,4) decrease quite fast with 
respect to variables x and y at infinity. Each bounded solution of the system (1.1) 
has the asymptotic 



tp(x,y) = F x a-(y) + o(l), y — ► -oo , 
ip{x,t) = F x a + (t) +o(l), y — > +oo . 

where F x denotes the diagonal matrix shift operator, such that for a vector function 
a (y) = col [on (y) , a 2 (y) , a 3 (y)] , F x a (y) = col [on (t + x) , a 2 (t + x) , a 3 (t - x)] 
The scattering operator S is defined by the equality 



(1.2) a + (y) = Sa_ (y) . 

This operator has the inverse 5 -1 and S = I + F and S^ 1 = I + G, where F = 
[Fij] i j =1 and G = [Gij\ i . =1 are matrix integral operators. To regenerate the system 
(1.1), i.e., the coefficients qi( i = 1, . . . , 4) it is sufficient to know Fi 3 , F 2 3, G31 and 
G32. The collection {-F13, -F23, G31, G32} is called the scattering data for the system 
(1.1) in the plane. The solution of the inverse problem reduces to the integral 
equation 
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/V I" (-+OC 
/ A (x, y, z)G(z — x, s + x)dz 
-co \-Jy 



B{x,y,r) 



where F 



+00 

y 

Pi 3 
P23 



B (x, y,z)F (z + x,s — x) dz 



and G = [ G 31 G 32 ] 



F(s + x,r — x)ds = F(y + x,T — x), 

(1.3) 

G(s + x,r — x)ds = G(y — x,T + x), 



The coefficients <fe( i — 1, ... ,4) of the system (1.1) are expressed in terms of 
the solution of (1.3) by means of the equations 



9i 

92 



(x, y) = -2A (x, y, y) , [ q 3 q 4 ] (x,y) = -2B (x,y,y) . 



3. Nonlinear evolution equation related to nonstationary 

DlRAC-TYPE SYSTEMS 

Let M and A be first order matrix operators in (0.2): 

M = a— + Q, A = t— — h P. 

ox ox 



Then the equation (0.2) becomes to the form 
(2.1) 



— - G— - Q — - T— - P 

dy dx ' dt dx 



0. 



Here r and P are third order square matrices. Let the matrix r be real and 
diagonal: r = diag{b\,b 2 ,b 3 ), b\ > b 2 > b 3 and the matrices Q and P obey the 
relation: [cr,P] = [t,Q]. 

wis 




Take P 



^r^9i 

^93 ^94 

solutions of the equations 







where the functions V12 and V21 are 



d 



d 



h - b 2 

-^-vi2 - -^-vi2 = — 9i94, 

oy ox 2 

d 



b 2 - bi 

_ W21 __ U21 = _ g2g3 , 

respectively. 

Then, the conditions (2.1) are reduced to the system of equations 
8 , d , 8 

7^91 + fc l^-<?l + k 2^~qi = "1292, 

at oy ox 

d , d , d 

^792 + ^1^-92+^2-^-92 = "2191, 

at oy ox 



d d d 

77793 + fa 77-93 + fa — q 3 = -W2194, 

at oy ox 

d , d , d 
TrM + fa^-qi + fa—qi = -v 12 q 3 , 
at oy ox 



(2.2) 
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where ki = -^2^, h = -^,fc 3 = - h -^M = -^f 3 -- 

Let us denote P = Jj — tJ^ — P. The Lax form (2.1) of the system of equations 

(2.2) enables us to apply to 1ST method for the integration. 

Lemma 1. Lettp be a solution of the Dirac-type system (1.1), whose the coefficients 
qi,i = 1, . . . , 4 satisfy Eqs. (2.2). Then the function ip = Pip also satisfy the system, 
(1.1). 

Proof. Let us apply the operator equation to ip: 

(LP - PL)V> = L(P^) - P(Lt/>) = 0. 
Taking into account that Li/> = 0, then 

L(PV) = 

follows from the last equation. It means that Pi[> is solution of system (1.1). □ 

The next theorem is true for the evolution of the scattering data. 

Lemma 2. Let the coefficients qi,i = 1, ... ,4 of the Dirac-type system (1.1) depend 
on t as a parameter and satisfy the system of equation (2.2). Besides that 

(2.3) vi 2 (x, ±oo ) = 0, V2i(x, ±oo) = 0. 

Then the kernels Fi 3 (y, r; t), F 23 (y, t; t), G 3 \(y, r; t), G 32 (y, t; t) of the integral 
operators F\ 3 , F23, G31, G32 corresponding to the scattering operator S for the sys- 
tem (1.1) in the plane satisfy the system of equations 

^-F 13 (y,T;t)~b 1 ^-F 13 (y,T;t)+b 3 -^-F 13 (y,T;t) = 0, 



(24) 



and 



d d d 

gl F 23 {V,T-,t) - b 2 —F 23 (y,r;t) + b 3 —F 23 (y,r;t) 



d d d 

—G 31 (y,T;t) + b 3 —G 31 (y,T;t)-b 1 —G 31 (y,T;t) = 0, 



d d d 

Q-G 32 {y,T;t)+b 3 —G 32 {y,T;t)-b 2 —G 32 {y,T;t) = 0. 



(2.5) 



Proof. By virtue of definition of the scattering operator S, from Lemma 1 we get 

(2.6) ip + — S(p_ 

where ip± = Po±, P = ^ - ra-§^. 
Since a + = Sa_ (see (1.2)), the equality 

(2.7) PS = SP 

follows from (2.6). 
Analogously, 

(2.8) PS _1 =S _1 P. 

From the matrix operator equation (2.7) it follows that the kernels of the integral 
operators F13 and F 23 satisfy system of equation (2.4). The similarly system of 
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equation (2.5) for the kernels of the integral operators G31 and G32 follows from 
the matrix operator equation (2.8). □ 

According to Chapter 2 and Lemma 2, let us give a procedure for the solution 
of the system (2.2) by 1ST method. 

Theorem 1. The system of equations (2.2) admits integration in the class of de- 
creasing functions that satisfy the condition (2.3) by the 1ST method. The solution 
qi,i = 1, . . . ,4 of the system (2.2) is determined by formulae 



<7i 

52 



(x, y; t) = -2A(x, y, y; t), [ q 3 q A ] (x, y; t) = -2B(x, y, y; t), 

where the vector functions A and B are the solutions of the integral equations 

A(x,y,r;t) - /"^ Jy°° A (x, y, z; t) G (z - x, s + x; t) dz F (s + x,t — x;t) ds 

= F(y + x,T-x;t), 

B (x,y,r;t) - Jy°° /"^ B (x,y, z;t) F (z + x,s - x;t) dz G (s + x,t — x;t) ds 

= G(y-x,T + x;t), 

and G = [ G31 G32 ] are satisfied the evolution equations 



where F — 



F13 

F23 

(24) and (2.5). 



Thus the system of nonlinear evolution equation (2.2) admits the integration 
by 1ST method. Let us denote by II the operator transforming the coefficients 
q(t) = col(qi(x,y;t),q2(x,y;t),q3(x,y;t),q4(x,y;t)) of the equation (1.1) to scat- 
tering data T(t) = col(-Fi 3 (j/, r; t), F 23 (y, r; t), G31O/, r; t), G 32 (y, r; t)) as follows: 

U:q(t)^T(t). 

Then the solution of the system (2.2) can be represent as 



(2.9) 



q{t) = U^e^T 



where T = T(0) i.e. T - co\(F 13 (y, r), F 23 (y, r), G 3 i(y, r), G 32 (y, r)) and 
bi£-h£ 



M 



b 2 ^-b 3 ^ 









The formula (2.9) is a specified form of solution q{t) of (2.2). 






b 2 4- - 63 7I7 
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